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For two-dimensional (2D) atomic Fermi gases in harmonic traps, the SO(2, 1) symmetry is broken 
by the interatomic interaction explicitly via the contact correlation operator. Consequently the 
frequency of the breathing mode lob of the 2D Fermi gas can be different from 2loq, with wo the 
trapping frequency of harmonic potentials. At zero temperature, we use the sum rules of density 
correlation functions to yield upper bounds for lob- We further calculate ujb through the Euler 
equations in the hydrodynamic regime. The obtained value of lob satisfies the upper bounds and 
shows deviation from 2loo which can be as large as about 8%. 
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Collective oscillation modes of atomic gases confined 
spatially in harmonic traps convey crucial information of 
the nature of the systems. The measurements of the fre- 
quency of the breathing mode in three dimensions have 
confirmed the superfluid hydrodynamics of Bose- Einstein 
condensates [T] and the universality of unitary Fermi 
gases [SJ. Moreover, it is pointed out in Ref. [3] that 
if the interatomic interaction satisfies scale invariance, 
the atomic gases possess a "hidden" SO(2, 1) symme- 
try. Such a symmetry dictates that the frequency of the 
breathing mode of the gases confined in a spherical har- 
monic potential with trapping frequency ioq must be 2loq. 

With the advent of experimental realization of two- 
dimensional atomic gases both of bosons [IHI] and 
fermions [8TU2|. Refs. [12] and [2] put into the context 
of quantum anomaly the point mentioned in Ref. [5] that 
if a contact pseudo-potential with a bare coupling con- 
stant g is used to model the short-ranged interatomic 
interaction, the two-dimensional atomic gases have the 
50(2, 1) symmetry at the classical level since g is dimcn- 
sionless. However, necessary renormalization of g intro- 
duces the scattering length a2D in two dimensions as a 
new low energy observable which characterizes the in- 
teratomic interaction. The classical 5*0(2, 1) symmetry 
ceases to hold at the quantum level due to the fact that 
a 2 D carries dimension. Consequently, the frequency of 
the breathing mode of the two-dimensional gases can 
be different from two times the harmonic frequency ujq 
of the spherical traps. On the other hand, recent ex- 
periment of two-dimensional Fermi gases carried out at 
temperature T rs 0.37Tp, with Tp the Fermi tempera- 
ture, did not observe substantial deviation of lob from 
2w at all [T5] . 

In this paper, we study the breathing mode of a two- 
dimensional (2D) Fermi gas of equal number of two 
species of fermionic atoms at zero temperature. Based 
on the short-rangedness of the interatomic interaction 
and the correlation structure at short distances in dilute 
Fermi gases, we give an alternative derivation to show 
that the interaction violates the SO(2, 1) symmetry ex- 
plicitly via the contact correlation operator, compared to 
the quantum field approach in Ref. [13]. We derive the 



sum rules of density correlation functions and use them 
to yield upper bounds for the frequency of the breathing 
mode ujb- We further calculate ub by the Euler equa- 
tions in the hydrodynamic regime and find that the ob- 
tained ujb satisfies the upper bounds and can be about 
8% bigger than 2wo in the unitary regime as shown in 
Fig. ([I]). Since in the high temperature limit atoms are 
nonintcracting and the system obeys the SO(2, 1) sym- 
metry trivially, ljb = 2cj - Our results at zero tempera- 
ture determine the typical order of magnitude how large 
ujb can deviate from 2wo at finite temperatures. Our 
results have the prospect to be verified by future experi- 
ment deep in the degenerate regime. 

Basic formalism. — Experiment produces two- 
dimensional Fermi gases by confining 6 Li atoms of two 
hyperfine states in a three-dimensional harmonic trap 
with a large trapping frequency ui z in the z direction 
[BTTTTj. Under the circumstances that all other relevant 
energy scales are much smaller than oj z , the motion 
of the atoms in the z direction is frozen; low energy 
dynamics occurs only in the xy plane. The system is 
described by the Hamiltonian 

H =H + Hhoi 

H ho =-mu>l E f ^> (!) 

with the operators f = {x,y}, p = {p X ifiy}i arL d m the 
atom mass and ujq the spherical trapping frequency in 
the xy plane. We denote the two species of the fermions 
by j" and J,. The attractive interaction potential U(r) 
has a short-range ro and gives rise to a shallow bound 
state with binding energy Ef, = —l/ma\ D , with the 2D 
scattering length a-iu r . (We take H = 1 through- 
out.) Since we are only interested in the long wave- 
length physics, the intra-species interactions have been 
neglected due to the Pauli exclusion principle. 

Experimentally the breathing mode can be excited by 
modulating the trapping frequency loq |15j . At zero tern- 



2 



perature, the breathing mode should bring about a sharp 
peak in the spectrum function of the density correlation 



fermions at short distances, obeys the adiabatic relation 

m 



x » = nJ2 \(f\Sd\g)\ 2 S(u - E, + E g ), (2) 
/ 

where \g) is the ground state, and O = J2i CT =t i anc ^ 
50 = O — (O). We use (...) to denote the expectation 
value over \g) . 

The structure of x"( w ) can be analyzed by calculating 
the equation of motion of the operator 0: 



The scale dilation generator is D = a (Pi,cr ' *i 



(3) 



Pi, CT )/2 since e^f^e 



,-i-yD _ „7 
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We continue 



to calculate the equation of motion for D: 



i~D(t)=[D,H] 



-2iH - AiH h 



d 2 r 



p(r) 
(4) 



Here the two-particle correlation function is p(r) — 
J2i j <K r — ^-t,t — If ^ ne interatomic interaction is 

scale invariant, i.e., U(e 1 r) = e~ 27 [/(r), the last term in 
Eq. Q is zero; given [D, O] = —2iO, the operators O, 
D and H form a closed algebra, which corresponds to 
a 50(2,1) symmetry [3]. In this case, the operator 
satisfies 



if 

df- 



o 



4 - 9 . 
— H - 4^0, 
m 



(5) 



from which one can read off that nonzero matrix elements 
(f\0(t)\g) oscillate with a frequency 2uj and conclude 
that x"(w) is a delta function centering at the frequency 
of the breathing mode lub — 2ujq. 

However, the real interatomic interaction is not scale 
invariant. To evaluate the last term in Eq. Q which 
breaks the 5*0(2, 1) symmetry, we note that due to the 
diluteness of atomic Fermi gases the two-particle corre- 
lation function has the asymptotic form [16 



p(r) = O0 2 (r) 



(6) 



for r -C d the mean interparticle spacing. The wave func- 
tion 4>{ r ) satisfies the two-body Schrodingier equation in 
the relative coordinates 



1 d d 

—r- — h mil (r) 

r dr dr 



(r) = 



(7) 



and is normalized such that <fr(r) = \n(r/a2D) f° r r ^ r o 
the range of U(r). The contact correlation operator 
C, which quantifies the correlation strength between 



(C) 



m d(H) 
2it dhia2D 



(8) 



at zero temperature. Similar to the manipulations em- 
ployed in Ref. [T7], we integrate the second term of the 
integrand in Eq. Q by parts, and by Eq. ^ find 



d 2 r 



n rr, N 9U(r 



dr 



cj) 2 (r) = -2ir/m. (9) 



Thus, we have 

[D,H] 



2iH - AiH ho + 2mC/m, 



(10) 



which agrees with the result derived by a quantum field 
approach in Ref. [M] . 

Sum rules. — With the contact correlation operator C 
breaking the 50(2, 1) symmetry explicitly, the breath- 
ing mode frequency ujb can be different from 2wo- To 
constrain the value of ujb, we use the sum rules, sg — 
I-oo dujx" ( UJ ) uji f 71 : to define the frequencies uigj-2 = 
\J sg/se-2', at zero temperature uig^-2 are upper bounds 
for ujb |18j . Specifically we consider 



S3 = ({{D,H],D}}, 
m z 

Sl =h[[sd,H],sd]) = -(d), 

2 m 



S-l 



1 d 

m duiQ 



(6). 



(ii) 

(12) 
(13) 



The third sum rule S3 can be evaluated as 
[[D,H],D] 

Ml -C I d 2 r Uu(r) - r^-r 9 ^ 

or or 



2 (r). (14) 



To proceed further, we assume a square well model po- 
tential U(r) = — Vo8(r Q — r) with Vq > and obtain 



d 2 r 



4U(r) 



d dU(r) 

r — r 

dr dr 



2 (r) 



47T 



m 



47rVbr 2 ln(r /a 2 D). 



(15) 



On the other hand, since the parameters Vq and ro are 
required to reproduce the low energy physical quantities, 
e.g., the scattering length a2D, in the limit ro/ci2D 0, 
we have 



ln(r /a 2 D) = - 



mV Q r 2 



From Eqs. (|14|)-(|16|) we obtain 

= — (H) 



3tt 



(C) 



(16) 



(17) 
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which also implies [D, C] — AiC where the factor 4 is 
the dimension of C/Q with f2 the system volume in two 
dimensions. 

By the virial theorem 



(H) 



muJ o(0)--(C), 



which can be derived from 



d 



1 



d 



dujQ 2 d In <i2D 
we cast W3 1 into the form 



- 1 )(H) = 0, 



U31 = 2uj q 



1 



2ir(C) 

W) 



,2,,,2 



1/2 



(18) 



(19) 



(20) 



On the other hand, using dimensional analysis, we find 
within the local density approximation the upper bound 



,_i = 2w /\/l - A 



with 



A = 



7r(2 + c)/91na 2D )(C') 
2m 2 a;2(0) 



(21) 



(22) 



We evaluate w 31 and u>i _i for the 2D Fermi gas 
trapped in the harmonic potential with the local den- 
sity approximation and using the interpolation of the 
equation of state obtained by the Monte Carlo simula- 
tion [19] as used in Ref. [12]. Figure shows the rela- 
tive deviations of the upper bounds W31 and ui\ 1 from 
2luq verse \n(kpa2D), where k 2 F — 2nnQ and uq is the 
total fermion density at the trap center. In the BEC 
limit a,2D + , fermions pair to form tightly bound 
bosonic molecules and these molecules are weakly in- 
teracting with each other. The ground state energy is 
(H) « — N/2ma\ D with N the total number of fermions. 
From Eq. pj), (C) N/2naj D , while (6) - N/^/mu^. 
The upper bound u^i diverges as 0,2 d — * + . The phys- 
ical reason of this divergence is that when the trapping 
frequency ujq is modulated, the operator O acting on \g) 
can disassociate the molecules. The integral of tu 3 x"(u) 
is dominated by the part at frequencies uj > \jma\ D 
(cf. Eq. (j2j). In the BCS limit a 2 D -> +00, the inter- 
action energy of the ground state is ~ — l/ln(<Z2D), and 
thus C ~ 1/ In 2 (122,0), while (O) is basically the value for 
nonintcracting fermions which is finite; W3.1 approaches 
2uj from above. 

The other upper bound does not suffer diver- 

gence in the BEC limit, which mathematically is due to 
lower powers of ui in Si than in S3. On the ground of 
dimensional analysis, (O) = /(mwodfj}, N)/lvq. We ex- 
press uji _i in terms of the function / as 



(23) 
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FIG. 1: (Color online) Relative deviations from 2loq versus 
ln(feFtS2D) for the upper bounds and W3,i, and lob cal- 

culated by the hydrodynamic approach . T he red curve is for 
8uj — — 2cjo evaluated from Eq. (20 1. The black curve 

in the inset is for Slo = L03.1 — 2loq from Eq. ( |21[ ). The blue 
squares are for <5cj = 0;^ — 2u>o evaluated from Eq. ( |30[ | for 
-2 < ln(fc F a2D) < 2. 



Since the cloud size increases as o,2d increases |12j . 
df(^,N)/d£_ > 0; is always bigger than 2wo, which 

implies A > 0. From the equation of state of homoge- 
neous gases obtained in Ref. [T5], one can deduce that 
Zdf(Z,N)/d£ becomes zero in the BEC and BCS limits 
and reaches a maximum at £ ~ 1. Correspondingly, as 
shown in Fig. ([I]), the bound w^-i approaches 2u>o in 
the BEC and BCS limits and shows a maximum in the 
unitary regime kp02D ~ E We find w^-i < 0*3,1 for any 
ln(fci?a2_o); the upper bound o^-i is more strict. 

Hydrodynamic equations. — The above sum rule results 
suggest that the frequency of the breathing mode ujb can 
differ from 2o>o by a significant amount in the unitary 
regime. The experiment [T5] showed that the degener- 
ate Fermi gas is in the hydrodynamic regime around the 
unitary limit fc^ 020 = 1< The small decay rate of the 
breathing mode measured there, which is primarily due 
to the anharmonicity of the trapping potential |20j . in- 
dicates that dissipation is negligible. We use the Eulcr 
equations 



On 
di 



+ V • (nv) = 0, 



■ v ■ Vv 



(24) 

-Vfi-VVho, (25) 



to calculate lob in the hydrodynamic regime at zero tem- 
perature. Here n and v are the density and velocity 
fields of the fermions, and /i(r) is the local chemical 
potential of the fermions. The harmonic potential is 



V ho (r) 



72. 



f=-row a§ D 



Previous applications of the Euler equations to atomic 
gases in three dimensions usually linearize the density 
n = n eq + Sn with n eq the equilibrium density distribu- 
tion, and treat both 5n and v as small quantities. If 
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one assumes Sn(r, t) 
equations give [21J 



e lut 5n(r), the linearized Euler 



-tu 2 Sn = V • 



» On'!*" 



(26) 



Here u eq is the equilibrium local chemical potential. It 
is tempting to solve Eq. ( 26 1 for the breathing mode of 
two-dimensional Fermi gases as an eigen-equation for Sn 
in the domain < r < Rtf with Rtf the Thomas- Fermi 
radius of the cloud at equilibrium where n eq (RTF) = 0. 
However, due to the attractive interatomic interaction, 
the part of the system close to the cloud edge is in 
the BEC regime. The equation of state for this part 
l/2ma| D ~ — n/\n(na 2 D ) 13 . Given that Sn for 



is fx 



the breathing mode should be spherically symmetric, the 
asymptotic form of Eq. (|26[) in this spatial region is 



dw 2 



_d_ 

dw 



1 



w In w 



Sn = 



(27) 



with w = 1 — r 2 /R TF . We discover Sn 



(-lnw) 1 / 2 

which diverges as w — > + [35]; no physical solution ex- 
ists. Mathematically speaking, the logarithm in the equa- 
tion of state of two-dimensional BECs makes r = Rtf 
an irregular singular point of Eq. ( 26 ) . We attribute the 



absence of physical solutions to Eq. (26) to the approx- 



imation of linearizing the Euler equations which must 
breaks down since the linearized solution Sn would be- 
come infinitely big close to the cloud edge. 

Physically the frequency of the breathing mode ojb is 
expected to be determined by the bulk of the cloud in- 
stead of the details close to the edge. To emphasize the 
property of the bulk, we adopt a variational approach 
and work with the Lagrangian 



L 



d 2 r 



1 



£ - nV h , 



On 



V • nv 



from which Eqs. (24 1 and (25) can be derived 



(28) 



The 



internal energy density £ is related to the local chemical 
potential u via d£/dn = u. The Lagrangian multiplier 

We 

; (r/A(t))/A(i) 2 



4> is used to ensure the continuity equation (|24|) 
assume a variational ansatz n(r,t) 



for the breathing mode which is exact if the interaction 
is scale invariant. Correspondingly Eq. (24) determines 
v(r,t) = (d\n\(t)/dt)r. Substituting the forms of n(r,t) 
and v(r, t) into Eq. (28), for small oscillations X(t) = 
1+S(t) with S(t) <C 1, we obtain the Lagrangian equation 



with 



d 2 S(t)/dt 2 +uJ 2 B S{t) = 0, 



4 J d 2 r n 2 eq (du eq /dn eq ) 



Jd 2 r 



(29) 



(30) 



Further manipulation yields 

u B = 2lo VT+~A, 



(31) 



from which one can tell u>b < -i for A > 



(cf. Eq. (21)). The physical meaning of Eq. (30) is man- 
ifest: its nominator is the average of the square of the 
sound velocity over the cloud and its denominator gives 
the square of the linear dimension of the cloud; lub is 
basically the rate how fast a density variation can propa- 
gate through the cloud with the averaged sound veloc- 
ity [23]. Note that the logarithmic singularity in the 
equation of state for two-dimensional BECs mentioned 
above is weak enough that the integral in the nominator 
of Eq. (|30|) converges. If the interaction is scale invariant, 



n eq (d^i eq /dn eq ) = fi eq ; given du eq /dr = -dV ho /dr, one 
retrieves lob = 2wo- 

With the local density approximation as used for u)\^-i 
and a^i, we evaluate u>b from Eq. (30) for —2 < 



hi(kpa2D) < 2 which is presumably in the hydrodynamic 
regime [15] . Figure ([T]) shows that the maximum devia- 
tion of ujb from 2w is about 8% at kpd2D ~ 1- The esti- 
mate of ujb in Ref. [14] . where a linearized Euler equation 
is solved with assuming a polytropic equation of state, is 
bigger than ours, and barely satisfies the upper bound 
0^1,-1 . Our calculation of ojb at zero temperature sets up 
the typical scale how large the difference between ujb a-ud 
2wo can be. It is natural to expect that when temperature 
rises up, ujb decreases from its value at zero temperature, 
and approaches 2ujq in the high temperature limit where 
the fermions are essentially noninteracting. Future ex- 
perimental advance into deep degenerate regimes has the 
prospect of suppressing the damping rate of the breath- 
ing mode, which is about 2% of 2w at T/T F sa 0.37 [IB] , 
and producing a definite detection of nonzero ujb — 2ujq. 

We acknowledge useful discussions with H. Zhai, 
Z. Chen, G. Baym, C.J. Pethick, M. Kohl, M. Rande- 
ria, E. Taylor, and M. Oshanii. We thank M. Oshanii for 
providing us with his notes on solving the Euler equa- 
tions for two dimensional gases. This work is supported 
in part by Tsinghua University Initiative Scientific Re- 
search Program, and NSFC under Grant No. 11104157. 



* Electronic address: huazhenyu2000@gmail.com 
[1] D.S. Jin, J.R. Ensher, M.R. Matthews, C.E. Wieman, 

and E.A. Cornell, Phys. Rev. Lett. 77, 420 (1996). 
[2] J. Kinast, S.L. Hemmer, M.E. Gehm, A. Turlapov, and 

J.E. Thomas, Phys. Rev. Lett. 92, 150402 (2004); 

M. Bartenstein, A. Altmeyer, S. Jochim, C. Chin, 

J.H. Denschlag, and R. Grimm, ibid. , 203201 (2004). 
[3] L.P. Pitaevskii and A. Rosch, Phys. Rev. A 55, R853 

(1997). 

[4] Z. Hadzibabic, P. Kriiger, M. Cheneau, B. Battelier, and 

J. Dalibard, Nature 441, 1118 (2006). 
[5] P. Kriiger, Z. Hadzibabic, and J. Dalibard, 

Phys. Rev. Lett. 99, 040402 (2007). 



5 



S. Tung, G. Lamporesi, D. Lobser, L. Xia, and E.A. Cor- 
nell, Phys. Rev. Lett. 105, 230408 (2010). 
X. Zhang, C. Hung, S. Tung, and C. Chin, Science 335, 
1070 (2012). 

K. Martiyanov, V. Makhalov, and A. Turlapov, 
Phys. Rev. Lett. 105, 030404 (2010). 
M. Feld, B. Frohlich, E. Vogt, M. Koschorreck, and M. 
Kohl, Nature 480, 75 (2011). 

B. Frohlich, M. Feld, E. Vogt, M. Koschorreck, W. Zw- 

erger, and M. Kohl, Phys. Rev. Lett. 106, 105301 (2011). 

A.T. Sommer, L.W. Cheuk, M.J.H. Ku, W.S. Bakr, and 

M.W. Zwierlein, Phys. Rev. Lett. 108, 045302 (2012). 

A. A. Orel, P. Dyke, M. Delehaye, C.J. Vale, and 

H. Hu, New J. Phys. 13, 113032 (2011). 

M. Olshanii, H. Perrin, and V. Lorent, 

Phys. Rev. Lett. 105, 095302 (2010). 

J. Hofmann, Phys. Rev. Lett. 108, 185303 (2012). 

E. Vogt, M. Feld, B. Frohlich, D. Pertot, M. Koschor- 
reck, and M. Kohl, Phys. Rev. Lett. 108, 070404 (2012). 

F. Werner and Y. Castin, |arXrv7l2 04.3204 

S. Zhang and A.J. Leggett, Phys. Rev. A. 79, 023601 
(2009). 



[18] L. Pitaevskii and S. Stringari, Bose- Einstein Condensa- 
tion, P. 168 (Oxford University Press, USA, 2003). 

[19] G. Bertaina and S. Giorgini, Phys. Rev. Lett. 106, 
110403 (2011). 

[20] S. Riedl, E.R.S. Guajardo, C. Kohstall, A. Altmeyer, 
M.J. Wright, J.H. Denschlag, R. Grimm, G.M. Bruun, 
H. Smith, Phys. Rev. A 78, 053609 (2008). 

[21] C.J. Pethick and H. Smith, Bose-Einstein Condensation 
in Dilute Gases, second edition (Cambridge University 
Press, Cambridge, 2008). 

[22] CM. Bender and S.A. Orszag, Advanced Mathematical 
Methods for Scientists and Engineers: Asymptotic Meth- 
ods and Perturbation Theory, P. 76 (Springer, New York, 
1999) 

[23] P.R. Zilsel, Phys. Rev. 79, 309 (1950); E. Taylor and 
A. Griffin, Phys. Rev. A, 72, 053630 (2005). 

[24] Since this work was largely completed, we became aware 
of E. Taylor and M. Randeria, arXiv: 1205. 1525 whose 
expr essi on for ojb at zero temperature is the same as 
Eq. j30l. 



